Zero-temperature phase diagram of Yukawa bosons 
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We study the zero-temperature phase diagram of bosons interacting via screened Coulomb 
(Yukawa) potential by means of the diffusion Monte Carlo method. The Yukawa potential is used 
cis a model interaction in the neutron matter, dusty plasmas and charged colloids. As shown by D. 
S. Petrov et al. [Phys. Rev. Lett. 99, 130407 (2007)], interactions between weakly bound molecules 
of heavy and light fermionic atoms are described by an effective Yukawa potential with a strength 
related to the heavy-light mass ratio M/m which might lead to crystallization in a two-dimensional 
geometry if the mass ratio of heavy-light fermions exceeds a certain critical value. In the present 
work we do a thorough study of the quantum three-dimensional Yukawa system. For strong in- 
teractions (equivalently, large mass ratios) the system experiences several phase transitions as the 
density is increased, passing from gas to solid and to gas phase again. Weakly interacting Yukawa 
particles do not crystallize at any density. We find the minimal interaction strength at which the 
crystallization happens. In terms of the two-component fermionic system, this strength corresponds 
to a heavy-light mass ratio of M/m ~ 180, so that it is impossible to realize the gas-crystal tran- 
sition in a conventional bulk system. For the Yukawa model of fermionic mixtures we also analyze 
the possibility of building molecular systems with very large effective mass ratios by confining the 
heavy component to a sufRciently deep optical lattice. We show how the effective mass of the heavy 
component can be made arbitrarily large by increasing the lattice depth, thus leading to a tunable 
effective mass ratio that can be used to realize a molecular superlattice. 

PACS numbers: 64.60.-i, 03.75.Ss, 71.10.Hf 



I. INTRODUCTION 

Recent advances in trapping and controlling ultracold 
dilute gases have permitted the realization of highly tun- 
able and extremely pure Fermi systems [1]. This has pro- 
vided new insight in the study of fundamental problems 
in condensed matter physics. For example, the original 
BCS theory[2] was developed to explain superconductiv- 
ity in metals, where the control over interactions and 
densities is very limited. However, in recent experiments 
with ultracold Fermi gases in the BCS-BEC crossover 
the strength of the interactions is controlled by external 
magnetic fields in the vicinity of a Feshbach resonance, 
while the geometry is tuned by means of magnetic or 
optical confinement. This has allowed, for instance, the 
measurement of the equation of state in the BCS-BEC 
crossover in high precision experiments [3]. Numerically, 
the best calculation of the zero-temperature equation of 
state is obtained in quantum Monte Carlo simulations [4] 

After the big success achieved with single species 
there is nowadays a growing interest in fermionic mix- 
tures. Quite recently, fermionic mixtures consisting of 
atoms with different masses have been realized experi- 
mentally [5, 6] and studied theoretically [7]. Novel phys- 
ical phenomena such as Efimov states [8], trimer and 
cluster formation might be observed [9] in these systems. 
The case of large mass imbalance is especially interest- 
ing, and mixtures of ^Li and ""'K are being investigated 
experimentally [5]. Even larger mass ratios are reached in 
mixtures of ^Li and "'^^^Yb[6]. In this article we present 
results for the phase diagram of Fermi mixtures as a func- 
tion of the mass ratio using quantum Monte Carlo meth- 



ods and determine how crystallization of this system can 
be realized. 

From the theoretical point of view, it was proposed 
in Ref. [10] that an effective Yukawa interaction, in- 
duced between heavy-light pairs of fermions, might lead 
to crystallization in quasi-two-dimensional systems. In 
this work we extend that discussion and analyze the pos- 
sibility of realizing a gas-crystal phase transition at zero 
temperature in three-dimensional (3D) systems. We ob- 
tain the phase diagram and discuss how large mass ratios 
have to be for reaching crystallization. 

The interest in the phase diagram of quantum Yukawa 
particles is rather old as the Yukawa potential has long 
been used, for instance, as a model for neutron mat- 
ter [11]. The Yukawa potential also describes interac- 
tions in dusty plasmas where charged dust particles are 
surrounded by plasma which introduces screening [12]. 
The Yukawa potential is often used as well as a model 
for suspensions of charged colloidal particles [13]. The 
classical finite temperature phase diagram has been ex- 
tensively studied [12, 13] while much less is known about 
the full quantum phase diagram. 

In the 70's, Ceperley and collaborators [14] used the 
diffusion Monte Carlo algorithm to estimate the zero- 
temperature phase diagram of the Yukawa Bose fluid. 
In their work the phase diagram was built assuming that 
the Lindemann ratio remains constant along the solid-gas 
coexistence curve, with the explicit value being evaluated 
only at a single point. In the present work we carry out 
a full study of the transition curve and present the phase 
diagram in terms of experimentally relevant densities and 
mass ratios of heavy to light fermions. The Lindemann 
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criterion prediction has turned out to be quite precise 
apart from the region of high densities. 



II. HAMILTONIAN 

Mixtures of fermions with different masses have been 
reaUzed recently in a new generation of experiments [5, 6] . 
The interactions can be tuned to ahow the formation 
of two-component molecules. The s-wave interactions 
within a single component are prohibited due to the Pauli 
principle. Yet. an ciffeetivc; interaction between same-spin 
fermions can be induced by the presence of the other 
component. The limit of large mass ratio has been an- 
alytically addressed in Ref. [10]. The effective interac- 
tion between heavy particles, which was obtained in the 
limit of large distances within first Born approximation, 
has the form of a screened Coulomb (Yukawa) potential. 
This leads to a description of the system in terms of a 
composite (molecular) bosonic gas interacting with an 
effective potential. The effective p-wave interaction be- 
tween heavy particles in BEC-BCS crossover, derived in 
Ref. [15], is no more of Yukawa form but rather present 
oscillations that increase in amplitude when going from 
BCS to BEC regime. In the present article we limit our- 
selves to considering the s-symmetry effective Yukawa 
interaction most relevant on the BEC side of the BEC- 
BCS crossover. 

We study a system of heavy fermions of mass M inter- 
acting among themselves and moving on a background 
of light fermions of mass m. The net effect induced by 
the movement of the light fermions can be characterized 
by a Yukawa potential, leading to the following effective 
Hamiltonian [10] describing the interaction between com- 
posite bosons formed by pairs of heavy and light atoms 

fr Y^A , 2fi^ exp(-2|r,-r,|/a) 



m 



where a is the atom-atom s-wave scattering length and 
Tj are positions of heavy atoms while the positions of 
light atoms have been integrated out. The ground-state 
properties of the system are then governed by two di- 
mensionless parameters, namely the gas parameter na^ 
and the mass ratio M/m. Equivalently, Hamiltonian (1) 
describes a bosonic system interacting via the screened 
Coulomb potential Vi„t(r) = gexp(— Ar)/r by mapping 
the charge to q = 2f^ /ma and the screening length to 
A-i =a/2. 

We calculate the ground-state properties correspond- 
ing to the Hamiltonian (1) by means of the diffu- 
sion Monte Carlo (DMC) algorithm [16]. This method 
solves stochastically the Schrodinger equation in imag- 
inary time providing the exact energy within control- 
lable statistical errors. The coexistence curves can then 
be traced by direct comparison of the energies of the 
solid and gas phases. The efficiency of the DMC method 
is greatly enhanced when importance sampling is used. 



This is done by multiplying the (unknown) ground-state 
wave function i/'(ri, . . . , r^r) by a guiding wave function 
■!/'T(ri, • • • 7 rw) and solving the equivalent Schrodinger 
equation for the product. As a result, the points in phase 
space where the guiding function is large get sampled 
more frequently and this improves convergence to the 
ground state. 

The properties of the gas phase are studied by con- 
structing the guiding function in a Bijl-Jastrow two-body 
product form VT(ri, • • • ,rjv) = Yli^j /2(|ri - r^l). We 
determine the optimal two-body Jastrow term /2(r) by 
solving the corresponding Euler Lagrange hypernettcd- 
chain equations[17] (HNC/EL) discarding the contribu- 
tion of the elementary diagrams. In this scheme the static 
structure factor S{k) that minimizes the variational en- 
ergy in the subspace of Jastrow wave functions has the 
form 



S{k) = 



(2) 



with t{k) = h^k"^ /2m and Vph{k) the so-called particle- 
hole interaction. Its Fourier transform FT\Vph{k)] = 
Vph{r) satisfies the following equation in coordinate space 

?l2 



Vph{r) = 9{r)V{r)+-\V^)\''+{9{ryi)uJi{r) , (3) 

where V{r) and g{r) are the bare two-body potential and 
the pair distribution function (the Fourier transform of 
S{k)), respectively. Finally, in momentum space the in- 
duced interaction 0Ji{k) becomes 

„K.) = -l.wMfi±ilM. (4) 

In this way, Eqs. (2), (3) and (4) form a set of nonlin- 
ear coupled equations that have to be solved iteratively. 
The Fourier transform of the resulting S{k) provides g{r) 
and, in this scheme, the optimal two-body Jastrow factor 
results from the corresponding HNC/0 equation 



f2{r) = Vg{r)(^~ 



N{r)/1 



(5) 



where N{r) is the sum of nodal diagrams, related to S{k) 
in momentum space by the expression N{k) = [S{k) — 

iY/s{k). 

The resulting wave function captures basic ingredi- 
ents coming both from the two- and many-body physics 
of the problem. On the other hand, the energy of 
the solid phase is obtained by using a Nosanow- Jastrow 

guiding wave function iI>t{yi, . . . ,vn) = nili/i(i"i " 
j.iatt.^ j^^jj.. _ 1^ ^ffli\l Gaussian one-body terms 

/i(ri — r-'^" ) = exp(— a(r — r'**" )^) describing the local- 
ization of particles close to the lattice sites r'^**'. The 
parameter a controls the localization strength and is op- 
timized by minimizing the variational energy. 

In order to find the energy in the thermodynamic limit, 
we carry out simulations of a system of N particles in 
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FIG. 1. (Color online) An example of finite-size dependence 
of the energy in the gas with periodic boundary conditions in 
truncated octahedron for M/m = 187 at two different den- 
sities na^ — 1.6 (upper set of data points) and na'^ — 0.192 
(lower set of data points). Symbols, DMC energy; lines, lin- 
ear fit to energy for large system sizes. Energies are scaled 
with the thermodynamic value Eextr, obtained in 1/A'^ — 
extrapolation. 



a box with periodic boundary conditions, and take the 
Umit N ^ oo while keeping the density fixed. In the sim- 
ulation of the crystal the number of particles should be 
commensurate with the box which restricts the allowed 
number of particles. For FCC packing the simulation 
box supports iV = 4z3 = 4, 32, 108, 256, .... In order to 
add more values we also use periodic boundary condi- 
tions on a truncated octahedron, which allows simula- 
tions with N = 2i^ = 2, 16, 54, 128, 250, 432, . . . particles 
with a larger effective volume of the simulation box and 
reduced anisotropy effects. Finally, the convergence is 
further improved by the Ewald summation technique[18] 
which we use in the calculations at large densities. 

In Fig. 1 we show two characteristic examples of the 
finite-size dependence of the energy at two different den- 
sities. For large enough system sizes, the energy is well 
fitted by a linear dependence in 1/A'^. For small number 
of particles the behavior is no longer linear, especially 
at large densities due to strong interparticle correlations. 
We find that system sizes of TV > 100 have to be used in 
order to ensure the linear regime at considered densities. 
The thermodynamic energy is then obtained as a result 
of a linear extrapolation 1/A^ — > 0. 



III. PHASE DIAGRAM 

An intrinsic property of Coulomb particles is to self 
assemble into a Wigner crystal at low densities and to re- 
main in a gaseous phase in the opposite limit, due to the 



long-range character of the interaction [19]. The Yukawa 
potential is similar to the Coulomb one at densities large 
enough for the interparticle distance to be much smaller 
than the screening length, which is fixed by the s-wave 
scattering length a. One then concludes that the Yukawa 
system stays in a gaseous phase at large densities. In the 
opposite regime of small densities, na'^ <^ 1, the interac- 
tion potential decays exponentially fast showing a short- 
range behavior that leads the system to a gaseous phase. 
For example, the FCC crystal of hard-sphere bosons of 
diameter Ug melts at density na^ ~ 0.24 [20]. The in- 
termediate regime na'^ « 1 is the most interesting one, 
as crystallization may or may not take place depending 
on the strength of the interaction, which in the current 
case of the Hamiltonian in Eq. (1) is governed by the 
mass ratio M /m. A relevant question then is what is 
the minimal mass ratio at which crystallization can be 
observed. 

In order to obtain an accurate description of the phase 
diagram, we study the finite size dependence and ex- 
trapolate the energy to the thermodynamic limit. The 
resulting energies of the gas and solid phases are then 
analyzed using the double-tangent Maxwell construction 
which provides the melting and freezing densities. The 
zero-temperature phase diagram parametrized in terms 
of the dimensionless density na^ and the mass ratio, is 
shown in Fig. 2. We find that for mass ratios smaller 
than the critical value M/m « 180 the gas phase is en- 
ergetically preferable at any density. On the other hand, 
for larger mass ratios there is always a gas-solid transi- 
tion at low densities and a solid-gas transition at large 
ones. Energetically, both the FCC and BCC lattices are 
possible in the solid phase. It is very difficult to discern 
numerically which packing is preferred as the energies in 
different crystalline phases are extremely close. Still, in 
the large potential energy limit, corresponding to a mass 
ratio M/m ^ 1, it is enough to compare the potential 
energy of the classical crystals with different packings. 
A simple, geometrical construction assuming that parti- 
cles are tightly tied to their equilibrium positions leads 
to a transition density na^ ~ 1.58. This prediction is 
depicted as a blue dashed line in Fig. 2. In the low- 
density limit we numerically find the value of the s-wave 
scattering length of the Yukawa potential (1) and fit 
it as as /a = 0.436 ln(Af/m) with accuracy below 1% in 
the region of interest. Note that a is the s-wave scat- 
tering length of fermionic particles which lead to the ef- 
fective bosonic Hamiltonian (1) while Os is the s-wave 
scattering length between bosonic Yukawa particles. For 
the sake of comparison we also plot in Fig. 2 the gas- 
solid transition line of hard spheres of size as given by 
M/m = exp(1.424/(ni/3o)). 

The figure also shows the results of Ceperley et al. [14] 
which were obtained by doing DMC calculations for three 
characteristic points in the phase diagram close to the 
solid-gas transition line. Overall, the agreement between 
that prediction and our results is good, the main differ- 
ences affecting the region of large density where Coulomb 
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FIG. 2. (Color online) Zero-temperature phase diagram of 
the Yukawa potential corresponding to the Hamiltonian in 
Eq. (1) in terms of the gas parameter na^ and the mass ratio 
M/m. Red symbols: transition point as obtained from the 
Maxwell double-tangent construction applied to the Monte 
Carlo data energies extrapolated to the thermodynamic limit; 
dashed line: critical density no? = 1.58 ... at which the en- 
ergy of perfect FCC and BCC packings are equal; dash-dotted 
line: prediction of Ceperley et al. [14] obtained by imposing a 
constant Lindemann ratio; short-dashed line: na^ — 0.24[20]. 



effects are strong. To our best knowledge this is the first 
time that the high-density quantum solid-gas phase tran- 
sition is observed in a simulation of Yukawa systems. 

In the case of the fermionic molecules, the resulting 
critical mass ratio is much larger than M/m « 13.6 for 
which the system is unstable due to formation of Efi- 
mov states [8]. The obtained phase diagram describes 
properties of metastable fermionic molecules while the 
true ground state corresponds to a many-body bound 
state. The stronger the effective interaction is (that is, 
the larger the mass ratio), the more distant are heavy 
fermions and the smaller the overlap with localized Efi- 
mov states is. 



IV. REACHING LARGE MASS RATIOS 

According to our results, the minimal mass ratio for 
which the crystalline phase can exist is M/m k, 180 and it 
is achieved at the somewhat large value of the gas param- 
eter na^ 0.3. At these densities the fermionic nature of 
the molecules becomes important as the Hamiltonian (1) 
is derived under the assumption that no? < 1/8[10]. Our 
bosonic model is expected to be reliable at smaller den- 
sities where the critical mass ratio is further increased. 

The mixtures of different fermionic atoms have already 
been successfully realized in experiments [6] but at sig- 
nificantly smaller mass ratios. Probably, the largest di- 



rectly achievable mass ratio currently is that of Yb and 
Li atoms, M/m — 29, which is still much smaller than 
the critical mass ratio needed to observe the formation 
of an ultracold crystal. 

An alternative way to realize a fermionic mixture with 
a large and variable mass ratio is to confine one of the 
components to an optical lattice. At low filling fraction 
the distances between atoms are large compared with the 
lattice spacing, and the separation of length scales allows 
the description of the movement of a particle in the lattice 
as that of a quasiparticle with an effective mass moving 
in a medium where the lattice is absent. In a deep lattice 
interactions between particles are much weaker than the 
confining energy and so, to a first approximation, one can 
consider that as the problem of a single particle diffusing 
in the lattice. 

An optical lattice created by counter-propagating laser 
beams imposes an external potential Matt, {x, y, z) — 
Vb (sin^ kx + bit? ky + sin^ kz) on every particle. The 
diffusion of a particle over a large distance is then gov- 
erned by the tunneling rate between neighboring sites. 
The diffusion is largely suppressed (and the effective 
mass greatly increased) when the amplitude of the op- 
tical lattice is large, i.e. when Vq ^ Er with Er = 
ti^k'^/2m the recoil energy. The excitation spectrum 
in the lowest band can be described by Bloch waves 
of quasi-momentum q and energy £o{q) = ffiwo — 
2 J (cos(j2:<i + cos g^d -I- cos (72C?) -I- . . . with d = ir/k the 
lattice constant [21]. At small momenta the spectrum is 
quadratic in q and can be interpreted as the spectrum 
eo{q) — Eq + fi^q'^ /2m* of a free quasiparticle with an 
effective mass to*. Within the lowest band approxima- 
tion the effective mass is inversely proportional to the 
hopping parameter J, 

TO* 1 Er 

— = —-r ■ (6) 

m TT^ J 

The tunneling is greatly suppressed in the deep opti- 
cal lattice limit Vq ^ Er- To better understand the 
contribution of the tunneling term in the present case, 
a semiclassical treatment within the Wentzel-Kramers- 
Brillouin(WKB) approximation can be used to calcu- 
late the tunneling probability. One finds that it is pro- 
portional to cx exp{— 2 J^^ dx y/2m{V{x) — E)/h}, 
where xi and X2 are the classical turning points. In the 
deep optical lattice limit the energy of the moving parti- 
cle is only slightly larger than the potential energy at a 
lattice site, and therefore V{r) — E V{r), so xi and X2 
correspond to the positions of two neighboring minima. 
The resulting integral can be easily evaluated and pre- 
dicts an exponential form J cx exp(— -y/Vb/^^r)- A more 
precise expression can be obtained from the width of the 
lowest band in the ID Mathieu-equation [21], yielding 



J 



zEr 



3/4 



exp ■ 



Vq_ 

Er 



1/2- 



(7) 



This expression, together with Eq. (6), provides an ana- 
lytic approximation for the effective mass m*. 



5 



In order to determine the dependence of m* on the 
lattice parameters in a non-perturbative way we evalu- 
ate the diffusion constant Z? of a real particle moving 
on the lattice and compare it to the diffusion constant 
Dq = /2m* of a free quasiparticle of effective mass 
m* . The diffusion constant is obtained by means of DMC 
propagation in imaginary time by measuring the mean- 
square displacement ((r(T) — r(0))^) = ((x(r) — a;(0))^) -I- 
{{y{T) - y{0)?) + {{zir) - ^(0))2) where r = {x, y, z) de- 
note particle coordinates. The diffusion constant is then 
extracted as I? = lim ?i((r(T) — r(0))'^)/(2Td), where 

r — >oo 

d = 3 is the system dimensionality. The resulting depen- 
dence of m* on the lattice amplitude is shown in Fig. 3. 
The figure shows the Monte Carlo prediction (solid line) 
compared with the approximation of Eq. (6) with J taken 
fromRef. [21] (circles) and from Eq. (7) (dashed line) . As 
it can be seen, there is an almost constant shift between 
m* obtained in the Monte Carlo simulation and Ref. [21] 
compared to Eqs. (6-7). We have found that the de- 
scription in the relevant region of interest is very much 
improved by subtracting a constant shift E'^^^ ~ —'ij^E^ 
from Vb in the argument of Eq. (7). This last prediction 
is shown by a thin line in Fig. 3 and provides a good 
approximation for Vq ^ lOE'r. 

One can understand these results in the following way: 
in the absence of the optical lattice the effective mass 
and the bare mass coincide, so m* = m. As the am- 
plitude Vo of the lattice is increased, the particle move- 
ment is slowed down and the effective mass increases. In 
the deep optical lattice limit the effective mass grows as 
m* /m (X exp(-\/Vo/-E'r) and so the ratio can be made 
arbitrarily large by increasing the amplitude Vq (for in- 
stance m* /m ^ 1000 at Vb/^-r = 40; see the inset in 
Fig. 3). This mechanism allows for increasing the mass 
of one of the two components while keeping the other one 
unaltered, so that the ratio M/m of the fermionic mix- 
ture can be made as large as desired when the mass of 
the heavy component is identified with the effective mass 
m*. Consequently, and according to the phase diagram 
shown in Fig. 2, there is a wide range of densities where 
one could find the system in the crystalline superlattice 
phase. Heights of optical lattices as large as (35 — 60)i?r 
are readily achieved in current experiments [22] and cor- 
respond to sufficiently large effective mass ratios for the 
crystallization to be realized. 

Both small density and large density transition lines 
are accessible for Yukawa interaction caused by screen- 
ing in dusty plasma, colloids and neutron matter. On 
the contrary, in two-component Fermi gas only the left 
part of the phase diagram can be realized since the ef- 
fective Yukawa interaction is valid only at low densities. 
In fact, the validity criterion for the interaction potential 
in Eq. (1) was studied in Ref. [10] and was found to be 
well satisfied for distances larger then r « 2a which leads 
to the condition pa^ ^1/8 when r is identified with the 
mean interparticle distance. In this way, for example, 
for pa^ = 0.1 and mass ratio M/m = 300 the system is 
expected to be in a crystalline form. Much larger effec- 




FIG. 3. (Color online) Effective mass as a function of the 
lattice amplitude Vo in units of the recoil energy E,-. Solid 
line: results obtained from the diffusion constant evaluated by 
propagation in imaginary-time; circles: lowest band approx- 
imation of Eq. (6) with values of J taken from [21]; dashed 
line: same results with J from the expansion in Eq. (7); dash- 
dotted line, same expansion with Vo shifted by —S/iEr- Inset: 
same results on a semi-logarithmic scale. 



five mass ratios can be achieved for realistic[22] lattice 
heights of (35 — 60)Er. We thus conclude that by using 
an optical lattice, a fermionic mixture of very different 
mass components can be used to test the phase diagram 
of the equivalent Yukawa model. 



V. SUMMARY AND CONCLUSIONS 

To summarize, in this work we have obtained the zero- 
temperature phase diagram of bosons interacting through 
Yukawa forces. We have used a diffusion Monte Carlo 
simulation starting from a very good approximation to 
the optimal variational ground-state wave function ob- 
tained by solving the corresponding Euler-Lagrange hy- 
pernetted chain equations. The resulting phase diagram 
is very similar to the one originally obtained by Ceper- 
ley and collaborators [14], although significant differences 
arise at large densities. The phase diagram shows that 
any fermionic mixture of pure elements will always be 
seen in gaseous form, as the mass ratios required for 
crystallization of weakly bound fermionic molecules are 
far beyond the ones that can be achieved in nature. Fi- 
nally, we investigate an alternative mechanism based on 
the confinement of one of the species to a deep optical 
lattice which exponentially increases its effective mass 
as a function of the confining amplitude. The resulting 
mass ratio of the mixture created in this way can then be 
tuned at will and could be used to check experimentally 
the predicted phase diagram both in the gas and crystal 
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(superlattice) phases. 



Grant No. 2009SGR-1003. G.E.A. acknowledges 
financial support from MEG (Spain). 



ACKNOWLEDGMENTS 



The work was sponsored by (Spain) Grant 
No. FIS2008-04403 and Generalitat de Gatalunya 



[1] F. Dalfovo, S. Giorgini, L. P. Pitaevskii, and S. Stringari 
Rev. Mod. Phys. 71, 463 (1999); S. Giorgini, L. P. 
Pitaevskii, and S. Stringari Rev. Mod. Phys. 80, 1215 

(2008). 

[2] J. Bardccn, L. N. Cooper, and J. R. Schrieffer Phys. Rev. 

108, 1175, (1957). 
[3] N. Navon, S. Nascimbene, F. Chevy, and C. Salomon 

Science 328, 729 (2010); M. J. H. Ku, A. T. Sommer, 

L. W. Cheuk, and M. W. Zwierlein, Science 335, 563 

(2012). 

[4] J. Carlson, S.-Y. Chang, V. R. Pandharipande, and K. 
E. Schmidt Phys. Rev. Lett. 91, 050401 (2003); S. Y. 
Chang, V. R. Pandharipande, J. Carlson, and K. E. 
Schmidt Phys. Rev. A 70, 043602 (2004); G. E. As- 
trakharchik, J. Boronat, J. Casulloras, and S. Giorgini 
Phys. Rev. Lett. 93, 200404 (2004); M. M. Forbes, S. 
Gandolfi, and A. Gezerhs Phys. Rev. Lett. 106, 235303 
(2011); J. Carlson, S. Gandolfi, K. E. Schmidt, and S. 
Zhang Phys. Rev. A 84, 061602(R) (2011). 

[5] M. Taglieber, A. C. Voigt, T. Aoki, T. W. Hansch, and 
K. Dieckmann Phys. Rev. Lett. ICQ, 010401 (2008); 

E. Wille, F. M. Spiegelhalder, G. Kerner, D. Naik, A. 
Trenkwalder, G. Hendl, F. Sclireck, R. Grimm, T. G. 
Tiecke, J. T. M. Wahaven, S. J. J. M. F. Kokkelmans, E. 
Tiesinga, and P. S. Julienne Phys. Rev. Lett. ICQ, 053201 
(2008); A. Voigt-C, M. Taglieber, L. Costa, T. Aoki, 
W. Wieser, T. W. Hansch, and K. Dieckmann Phys. 
Rev. Lett. 102, 020405 (2009); F. M. Spiegelhalder, A. 
Trenkwalder, D. Naik, G. Hendl, F. Schreck, R. Grimm 
Phys. Rev. Lett. 103, 223203 (2009); T. G. Tiecke, M. R. 
Goosen, A. Ludewig, S. D. Gensemer, S. Kraft, S. J. J. M. 

F. Kokkelmans, and J. T. M. Walraven Phys. Rev. Lett. 
104, 053202 (2010); L. Costa, J. Brachmann, A. C. Voigt, 

C. Hahn, M. Taglieber, T. W. Hansch, and K. Dieckmann 
Phys. Rev. Lett. 105, 123201 (2010); A. Trenkwalder, C. 
KohstaU, M. Zaccanti, D. Naik, A. I. Sidorov, F. Schreck, 
and R. Grimm Phys. Rev. Lett. 106, 115304 (2011); A. 
Ridinger, S. Chaudhuri, T. Salez, U. Eismann, D. R. Fer- 
nandes, K. MagaJhaes, D. Wilkowski, C. Salomon, and 
F. Chevy Eur. Phys. J. D 65, 223 (2011). 

[6] V. V. Ivanov, A. Khramov, A. H. Hansen, W. H. Dowd, 
F. Miinchow, A. O. Jamison, and S. Gupta Phys. Rev. 
Lett. 106, 153201 (2011); H. Hara, Y. Takasu, Y. Yar- 
maoka, J. M. Doyle, and Y. Takahashi Phys. Rev. Lett. 
106, 205304 (2011). 

[7] A. Gezerlis, S. Gandolfi, K. E. Schmidt, and J. Carlson, 
Phys. Rev. Lett. 103, 060403 (2009); I. Bausmerth, A. 
Recati, and S. Stringari, Phys. Rev. A 79, 043622 (2009); 

D. Blume and K. M. Daily Phys. Rev. Lett. 105, 170403 
(2010). 

[8] D. S. Petrov Phys. Rev. A 67, 010703(R) (2003); K. Hel- 
frich, H.-W. Hammer, and D. S. Petrov Phys. Rev. A 



81, 042715 (2010); Guan-Qiang Li and Ping Peng Phys. 
Rev. A 83, 043605 (2011); Tao Yin, Peng Zhang, and Wei 
Zhang Phys. Rev. A 84, 052727 (2011); Yujun Wang, W. 
Blake Laing, Javier von Stecher, and B. D. Esry Phys. 
Rev. Lett. 108, 073201 (2012). 
[9] T. Kraemer, M. Mark, P. Waldburger, J. G. Danzl, C. 
Chin, B. Engeser, A. D. Lange, K. Pilch, A. Jaakkola, 
H.-C. Nagerl, and R. Grimm Nature 440, 315 (2006); S. 
Knoop, F. Ferlaino, M. Mark, M. Berninger, H. Schobel, 
H.-C. Nagerl, and R. Grimm Nature Phys. 5, 227 (2009); 
M. Zaccanti, B. Deissler, C. D'Errico, M. Fattori, Jona- 
Lasinio M., S. Miiller, G. Roati, M. Inguscio, and G. 
Modugno Nature Phys. 5, 586 (2009); G. Barontini, C. 
Weber, F. Rabatti, J. Catani, G. Thalhammer, M. In- 
guscio, and F. Minardi Phys. Rev. Lett. 103, 043201 
(2009); J. R. Williams, E. L. Hazlett, J. H. Huckans, R. 
W. Stites, Y. Zhang, and K. M. O'Hara Phys. Rev. Lett. 
103, 130404 (2009); A. N. Wenz, T. Lompe, T. B. Otten- 
stein, F. Serwane, G. Ziirn, and S. Jochim Phys. Rev. A 
80, 040702 (2009); S. Knoop, F. Ferlaino, M. Berninger, 
M. Mark, H.-C. Nagerl, and R. Grimm Journal of Physics 
Conference Series 194, 012064 (2009). 
Shuta Nakajima, Munekazu Horikoshi, Takashi 
Mukaiyama, Pascal Naidon, and Masahito Ueda 
Phys. Rev. Lett. 105, 023201 (2010). 
[10] D. S. Petrov, G. E. Astrakharchik, D. J. Papoular, C. 
Salomon, and G. V. Shlyapnikov Phys. Rev. Lett. 99, 
130407 (2007). 

[11] G. Baym, and C. J. Pethick Ann. Rev. Nucl. Sci. 25, 27 
(1975); Hikaru Kawamura Prog. Theor. Phys. 66, 772 
(1981). 

[12] S. Hamaguchi, R. T. Farouki, and D. H. E. Dubin, 
J. Chem. Phys. 105, 7641 (1996); S. Hamaguchi, R. 
T. Farouki, and D.H.E. Dubin, Phys. Rev. E 56, 4671 
(1997); P. K. Shukla and K. Avinash Phys. Rev. Lett. 
107, 135002 (2011). 

[13] K. Kremer, M. O. Robbins, and G. S. Crest, Phys. Rev. 
Lett. 57, 2694 (1986); O. Robbins, K. Kremer, and G. S. 
Grest, J. Chem. Phys. 88, 3286 (1988); R. O. Rosenberg 
and D. Thirumalai, Phys. Rev. A 36, 5690 (1987); E. J. 
Meijer and D. Frenkel, J. Chem. Phys. 94, 2269 (1991); 
M. Heinen, P. Holmqvist, A. J. Banchio, and Gerhard 
Nagele J. Chem. Phys. 134, 044532 (2011); J. Gapinski, 
G. Nagele, and Adam Patkowski J. Chem. Phys. 136, 
024507 (2012). 

[14] D. M. Ceperley G. V. Chester, and M. H. Kalos Phys. 
Rev. D 13, 3208 (1976); D. M. Ceperley G. V. Chester, 
and M. H. Kalos Phys. Rev. B 17, 1070 (1978). 

[15] Y. Nishida, Phys. Rev. A 79, 013629 (2009). 

[16] For a general reference on the DMC method sec, e.g., J. 
Boronat and J. Casulleras Phys. Rev. B, 49, 8920 (1994). 



7 



[17] E. Krotschcck Microscopic Quantum Many-Body Theo- 
ries and Their Applications, edited by J. Navarro and A. 
Polls (Springer, Berlin, 1998). 

[18] O. N. Osychenko, G. E. Astrakharehik, and J. Boronat 
Molecular Physics 110, 227 (2012). 

[19] E. Wigner Phys. Rev. 46, 1002 (1934). 

[20] J. P. Hansen, D. Levesque, and D. Shiff Phys. Rev. A 
3, 776 (1971); M. H. Kalos, D. Levesque, and L. Vcrlet 
Phys. Rev. A 9, 2178 (1974); A. R. Denton, P. Nielaba, 
K. J. Runge, and N. W. Ashcroft Phys. Rev. Lett. 64, 
1529 (1990). 



[21] I. Bloch, ,J. Dalibard, and W. Zwerger Rev. Mod. Phys. 

80, 885 (2008). 

[22] K. Winkler, G. Thalhammer, F. Lang, R. Grimm, J. 
Hecker Denschlag, A. J. Daley, A. Kantian, H. P. Buchler, 
and P. Zoller Nature (London) 441, 853 (2006); G. Thai- 
hammer, K. Winkler, F. Lang, S. Schmid, R. Grimm, and 
J. Hecker-Denschlag Phys. Rev. Lett. 96, 050402 (2006); 
J. Catani, G. Lamporcsi, D. Naik, M. Gring, M. Inguscio, 
F. Minardi, A. Kantian, and T. Giamarchi Phys. Rev. A 
85, 023623 (2012). 



